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29ABSTRACT 

30Continuous measurements of sap flow have been widely used to measure water flux 
31through tree stems and branches. However, these measurements lack the resolution 
32necessary for determining fine-scale, leaf-level responses to environmental variables. 
33We used the heat ratio method to measure sap flow rates through leaf petioles and 
34leaflet petiolules of saplings of the tropical tree Tabebuia rosea (Bignoniaceae) to 
35determine how leaf and leaflet sap flow responds to variation in light and vapor pressure 
36deficit (VPD). We found that in the morning sap flow rates to east-facing leaves 
37increased 26 minutes before adjacent west-facing leaves. Although leaves had higher 
38integrated sap flow than their largest leaflet, this difference was not proportional to the 
39difference in leaf area, which could be due to lower conduit area in petiolules than in 
40petioles. In contrast to measurements on main stems, integrated daily sap flow was 
41 negatively correlated with daily mean VPD. Furthermore, leaves exhibited previously 
42undescribed patterns of hysteresis in the sap flow-VPD and sap flow-PAR relationships. 
43When hysteresis in the sap flow-PAR relationship was clockwise, the sap flow-VPD 
44relationship was also clockwise; however, when hysteresis in the sap flow-PAR 
45relationship was counterclockwise, the sap flow-VPD relationship displayed an 
46intersected loop. These pattern differences highlight how substantially leaf-level 
47processes may vary within a canopy and how leaf-level processes may not scale 
48predictably to the stem level. 
49 
50 
51 
52 
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53INTRODUCTION 

54 Approximately 90% of all water converted from the liquid phase to the vapor 
55phase in terrestrial ecosystems moves through plants (Jasechko et al. 2013); in the 
56tropics this amounts to an estimated 32x1 0 15 kg of water per year (Hetherington and 
57Woodward 2003). Almost all of this water transits through plant leaves. Understanding 
58how leaves respond to abiotic drivers is important for modeling efforts at scales from 
59leaves to landscapes (Jarvis and McNaughton 1986). At the leaf level, knowing which 
60drivers impact transpiration most under daily and seasonally varying conditions is critical 
61to understanding what may limit the distribution and abundance of species across the 
62globe. 

63 Various sap flow methods are commonly used to estimate almost continuously 
64tree responses to environmental conditions for extended time periods (Marshall 1958, 
65Granier 1985, Burgess et al. 2001, Vandegehuchte and Steppe 2012). These 
66measurements are most often made on boles or large branches of trees and can be 
67used to estimate canopy-level responses to changing environmental conditions (Oren et 
68al. 1999a, Ewers and Oren 2000, Traver et al. 2010). Despite recent technical 
69advancements in using sap flow measurements in stems to estimate canopy processes, 
70a number of problems remain. First, resistance and capacitance in the hydraulic 
71 pathway creates time lags in water movement between different points in the hydraulic 
72pathway (Kocher et al. 2013). For tropical trees, lags in sap flow between branches in 
73the canopy and the stem base can approach an hour (Meinzer et al. 2004). Second, 
74different parts of the canopy respond to environmental conditions largely independently 
75(Brooks et al. 2003), such that sap flow measurements on branches or boles may 
76provide only an average response of many leaves or branches. For example, east- 
77facing branches of Sequoiadendron giganteum reached their maximum daily sap flow 
78rates 6 hours before west-facing branches at the same height (Burgess and Dawson 
792008). Although sap flow measurements on boles and branches have provided useful 
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80estimates of whole-tree water use (Wullschleger et al. 1998), their utility for describing 
81 leaf-level processes can be limited by a variety of factors including time lags, 
82capacitance, hydraulic resistance, and variation in these factors along the root-to-leaf 
83continuum. 

84 Rarely have researchers attempted to measure sap flow rates through petioles of 
85individual leaves (Sheriff 1972). Recently, Clearwater (2009) adapted the heat ratio 
86method (HRM; Burgess et al. 2001) originally used for measuring sap flow through large 
87stems to measure sap flow through small diameter stems, fruit pedicels, and leaf 
88petioles. Slight variations of this method have proven useful in measuring sap flux 
89through petioles under field conditions in the neotropics (Roddy and Dawson 2012, 
902013, Goldsmith et al. 2013) and through stems of anatomically and phylogenetically 
91diverse species of the South African fynbos flora (Skelton et al. 2013). These studies 
92show that measuring sap flow directly adjacent to transpiring leaves can deepen our 
93understanding of how leaves respond to variation in environmental conditions across a 
94range of timescales. Placing sensors in close proximity to the transpiring leaves has the 
95advantage of fine-scale measurements akin to leaf gas exchange without the 
96disadvantage of enclosing leaves in a cuvette that removes the leaf boundary layer and 
97otherwise modifies the leaf microenvironment. 

98 Variation in sap flux is influenced by a variety of environmental conditions, 
99including soil water availability, vapor pressure deficit (VPD), and solar radiation, and 
lOOdiurnal patterns may also vary seasonally (e.g. O'Grady et al. 1999, 2008, Zeppel et al. 
1012004). Over diurnal cycles, a change in an environmental variable in the morning does 
102not always produce an equivalent response in sap flow as it does in the afternoon. 
103Such a pattern is termed hysteresis and has been commonly observed in the sap flow 
104responses to light and VPD. For example, at a given VPD that occurs both in the 
105morning and again in the afternoon, sap velocity is higher in the morning (when VPD is 
106increasing) than in the afternoon (when VPD is decreasing), creating a clockwise 
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107pattern of hysteresis throughout the day (Meinzer et al. 1997, O'Grady et al. 1999, 
108Zeppel et al. 2004) that is consistent with hysteresis in canopy gas exchange (Takagi et 
109al. 1998). Despite the ubiquity of hysteresis in sap flow data, its causes are rarely 
HOdiscussed. Hysteresis in a relationship indicates that factors other than the primary 
llldescriptor variable are constraining the response variable. For the relationship between 
112sap flow and VPD, it is thought that hysteresis results from variation in hydraulic 
113capacitance, resistance, or stomatal sensitivity to VPD (O'Grady et al. 1999). In the 
114simplest case, no hysteresis in the sap flow-VPD relationship would mean that an 
115increase and a decrease in VPD produce equivalent responses in sap flow and that 
116there are no other factors influencing the sap flow response to VPD. However, a 
117number of factors could cause deviation in the morning and the afternoon from this 
118scenario of no hysteresis. First, trees often supply their morning transpiration from 
119capacitive stores, which could elevate morning sap flow above that observed if there 
120were no capacitance (Cowan 1972, Goldstein et al. 1998, Meinzer et al. 2004, 2008). 
121Second, as stem water potential declines throughout the day, resistance in the hydraulic 
122pathway increases, as is commonly observed in vulnerability curves (Meinzer et al. 
1232009), which could depress afternoon sap flow below that if there were no resistance. 
124Both of these hydraulic factors, resistance and capacitance, could jointly be responsible 
125for causing hysteresis in the sap flow-VPD relationship. Morning transpiration may draw 
126largely on stored water, and as this hydraulic capacitor discharges, resistance may 
127become important in depressing afternoon sap flow. 

128 In contrast to VPD hysteresis, causes of hysteresis in the sap flow-light 
129relationship are less clear. While the sap flow-VPD relationship commonly exhibits a 
BOclockwise pattern of hysteresis, the sap flow-light relationship generally exhibits a 
131counterclockwise pattern (Meinzer et al. 1997, Zeppel et al. 2004). Zeppel et al. (2004) 
132argued that counterclockwise hysteresis in the sap flow-light relationship results from 
133the combination of (1) the difference in timing between peak light and peak VPD and (2) 
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134different stomatal responses to light and to VPD. Because VPD reaches its daily peak a 
135few hours after light reaches its daily maximum at solar noon, sap velocity will be higher 
136in the afternoon, when VPD is higher and stomata are fully open. Zeppel et al. (2004) 
137argue that stomatal conductance saturates at relatively low light levels in the morning, 
138and that above this saturating light level, VPD becomes the predominant driver of 
139transpiration and sap flow. These responses probably vary between leaves acclimated 
140to different microenvironments (e.g. between sun- and shade-leaves). 
141 Using sap flow measurements on main stems of canopy trees to test these 
142hypotheses for the causes of hysteresis are thus fraught with potential problems that 
143focusing on leaves may circumvent. For leaves, sap flow responses to VPD are often 
144similar to those for stems, although under some conditions leaves show different 
145patterns (Roddy and Dawson 2013). In addition to clockwise hysteresis in the 
146responses to VPD, sap flow through petioles sometimes exhibits an intersected loop (or 
147'figure-eight') pattern in response to diurnal variation in VPD. Determining the sap flow 
148responses to environmental variables of individual leaves provides an opportunity to 
149better elucidate important dynamics of plant water use. Furthermore, incorporating 
150explicit measurements of sap flow to individual leaves could help to improve upon 
151methods for scaling up to whole canopy processes. 

152 In the present study, we measured sap flow rates through petioles and petiolules 
153of saplings of the tropical tree Tabebuia rosea (Bignoniaceae) to understand how sap 
154flow responds to variation in light and VPD. Because the figure-eight pattern of 
155hysteresis in the sap flow-VPD relationship reported by Roddy and Dawson (2013) may 
156result from an interaction with light, we also measured photosynthetically active 
157radiation (PAR) levels on each leaf or leaflet to determine the conditions under which 
158different patterns of hysteresis may occur. We were particularly interested in examining 
159the differences in sap flow patterns between adjacent leaves and between leaves and 
160leaflets because different microenvironmental conditions may cause sap flow patterns to 
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161differ between leaves on the same stem. Furthermore, differences in sap flux through 
162petioles and petiolules may reflect variation in hydraulic architecture. If the hydraulic 
163pathway constricts downstream, then sap velocity must increase as it moves towards 
164leaflets. Our results highlight how measuring sap flow rates to individual leaves could 
165deepen our understanding of the linkages between hydraulic architecture and plant 
166water use. 
167 

168METHODS 
169Plant Material 

170 Tabebuia rosea (Bignoniaceae) grows to become a canopy tree in the lowland 
171forests of central Panama. While adults are often deciduous, seedlings are evergreen 
172with five palmate leaflets of varying size encircling the petiole. Plants were grown from 
173seed in 20-liter, insulated pots outdoors under a glass roof at the plant growth facilities 
174of the Smithsonian Tropical Research Institute in Gamboa, Panama, until a few days 
175before sap flow measurements were begun. When the tops of the plants were ~70 cm 
176above the soil surface, they were transferred to a glass chamber to protect them from 
177strong afternoon winds. The lower ~1 m of this chamber was made of cement painted 
178white, and doors on east- and west-facing sides of the chamber were left ajar to allow 
179air circulation. At the beginning of the sap flow measurements, the tops of the plants 
180were even with the top of the cement wall at the base of the chamber, and during the 
181course of the experiment two new sets of leaves were produced. Pots were kept well- 
182watered except for one week when water was withheld to determine how sap flow rates 
183would respond to declining soil water. This week coincided with a dramatic increase in 
184VPD. Sensors were installed when plants were approximately eight months old, a few 
185days after transferring them to the glass chamber. 
186 

187Sap flow measurements 
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188 On each measured leaf, sap flow sensors were installed on the leaf petiole and 
189on the petiolule of the middle, largest leaflet. On each plant, two adjacent leaves on 
190opposite sides of the plant were chosen for measurement. At the time of installation, 
191these leaves were the newest, fully-expanded leaves on each plant. Plants were 
192positioned so that the axis defined by the two measured leaves on each plant were 
193oriented east-west. Of the total 12 sensors installed, five failed, leaving two sensors on 
194petiolules and five sensors on petioles. 

195 Sap flow sensors and measurements were based on the design and theory of 
196Clearwater (2009) with some slight modifications described previously (Roddy and 
197Dawson 2012, 2013) and again briefly here. Sensors were constructed from a silicone 
198backing and were connected to 10 cm leads with Molex connectors that were then 
199connected by 10 m leads to an AM16/32 multiplexer and CR23X datalogger (Campbell 
200Scientific Inc., Logan, UT). Sensors were held in place with parafilm, and sensors and 
201connections were insulated with multiple layers of bubblewrap and aluminum foil at least 
2022 cm above and below the sensor. Consistent with previous applications of the HRM, 
203we measured initial temperatures for 10 seconds prior to firing a 4-second heat pulse, 
204monitored temperatures every 2 seconds for 200 seconds after the heat pulse, and 
205initiated the measurement routine every 10 minutes. 

206 The heat pulse velocity, v h (cm s 1 ), was calculated from the temperature ratio as: 

207 



h ~ 

X 



5T, 



\ v± 2 J 

208where v h is the heat pulse velocity in cm s 1 , k is the thermal diffusivity (cm 2 s 1 ), x is the 
209distance from the heater to each of the thermocouples (cm), and 5T 1 and 5T 2 are the 
210temperature rises (°C) above and below the heater, respectively (Marshall 1958, 

211Burgess et al. 2001 , Clearwater et al. 2009). We estimated the thermal diffusivity as: 

212 
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8 



8 



Downloaded from http://biorxiv.org/on September 18, 2014 

213where t m is the time (seconds) between the heat pulse and the maximum temperature 
214rise recorded xcm above or below the heater under conditions of zero sap flow 
215(Clearwater et al. 2009). We measured t m every morning before dawn when 
216atmospheric vapor pressures were lowest (between 0500 and 0600 hrs). At this time, 
217the vapor pressure deficit was almost always below 0.3 kPa, and therefore we assumed 
218v„ was approximately zero. Thermal diffusivity, k, was calculated for each thermocouple 
219(upstream and downstream) from these predawn measurements of t m , averaged for 
220each sensor, and used to calculate v h from the heat ratios for the subsequent 24 hours. 
221 Measurements of /con nights with VPD always above 0.3 kPa were discarded and 
222replaced with the most recently measured /cwhen VPD < 0.3 kPa. We estimated the 
223temperature ratio under zero-flow conditions by excising petioles and petiolules above 
224and below the sensor at predawn at the end of the experiment, greasing the cut ends, 
225placing the segments in a darkened box, and recording the temperature ratios for the 
226subsequent ~4 hours. The average of these zero-flow temperature ratios corresponded 
227very well with the temperature ratios recorded predawn under low VPD (less than ~0.3 
228kPa) conditions. The sensor-specific average temperature ratio under zero-flow 
229conditions was subtracted from all calculated heat ratios. This corrected heat ratio was 
230then used to calculate v h . 
231 

232Measurements of light and vapor pressure deficit 

233Light measurements were made using S1787 photodiodes (Hamamatsu Photonics, 
234Hamamatsu City, Japan). Photodiodes were connected to 15 cm long copper wires with 
235Molex connectors and then to 10 m leads, which were connected to a CR5000 
236datalogger measuring in differential mode. Circuits created by each photodiode were 
237closed with a 100 Ohm resistor. Photodiodes were installed just above each leaflet with 
238a sap flow sensor, and the photodiode was positioned to be parallel to the axis of the 
239central vein of the leaflet. Light measurements were made every minute and averaged 
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240every 10 minutes. Voltage measurements from the photodiodes were converted to PAR 
241based on a calibration of all photodiodes against a PAR sensor (LI-SB190, LiCor 
242Biosciences, Logan, UT), during which time all sensors were situated adjacent to each 
243other in a clearing that received full sunlight. 

244 Vapor pressure deficit was calculated from temperature and relative humidity 
245measurements made every 10 minutes with a HOBO U23 datalogger (Onset Computer 
246Corp., Bourne, MA) that was housed in a covered, white, PVC, Y-shaped tube and hung 
247level with the tops of the plants. 

248 

249Data analysis 

250 All analyses were performed using R (R Core Team 2012). Raw velocity 

251 measurements were processed following previously published methods (Roddy and 
252Dawson 2012, 2013, Skelton et al. 2013). Measurements of v h were smoothed using 
253the 'loess' function, which fits a polynomial to a subset of the data in a moving window 
254of 35 points. 

255 For analyses of structure (leaf vs. leaflet) or aspect (east- vs. west-facing leaves), 
256sap flow measurements from individual sensors in each group were averaged. To 
257estimate the total sap flow during the day and night, we integrated the time course of vh 
258measurements for each day and night using the 'auc' function in the package MESS, 
259which calculates the area under the curve using the trapezoid rule. Daytime was 
260defined as being between 600 and 1800 hours, which corresponded to morning and 
261evening twilight. To minimize the effects of nocturnal refilling, we defined nighttime as 
262being between 100 and 600 hours, which assumed that diurnal water potential declines 
263had mostly recovered within seven hours after sunset. To analyze the effects of VPD on 
264integrated sap flow rates, we linearly regressed integrated sap flow against mean VPD. 
265ln all regressions for leaves and leaflets in the day and in the night, the linear model 
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266was determined to be as good or better than both the logarithmic and power functions 
267by comparing the residual standard errors. 

268 Differences in the timing of morning sap flow between east- and west-facing 
269leaves were compared at a critical v h of 1 .5 cm hr 1 . We chose this critical value 
270because it was higher than any measured nighttime velocities and lower than most 
271daytime velocities. We estimated the time at which v h = 1 .5 cm hr 1 by assuming a linear 
272relationship (y = mx + b) between the two sequential morning measurements that 
273spanned v h of 1 .5 cm hr 1 . To compare east- versus west-facing leaves and leaves 
274versus leaflets, we used linear mixed effect models with day as the random variable, 
275which accounts for repeated measures. 
276 

277RESULTS 

278 Daily maximum VPD varied from 2.8 kPa to 6.9 kPa during the experiment, while 
279daily maximum PAR at the top of the canopy varied from 1125 to 1640 ji/mol nr 2 sec 1 . 
280The daily maximum sap flow rate through petioles varied between 1 .4 cm hr 1 to 4.5 cm 
281hr 1 . This lowest daily maximum v h occurred at the end of a week without water, during 
282which time five of the seven hottest, driest days occurred. Nighttime v h through petioles 
283varied throughout the study, but was always below 1 .0 cm hr 1 and below 0.5 cm hr 1 on 
284all but seven nights. Overall, thermal diffusivity, k, ranged from 0.00136 to 0.00170 cm 2 
285s 1 . There were slight differences in k between sensors, but /cwas relatively constant 
286throughout the experiment for each sensor, consistent with previously reported values 
287for k from a diverse set of plant structures and species (Clearwater et al. 2009, Roddy 
288and Dawson 2012, 2013, Skelton et al. 2013). 

289 On every morning, sap flow rates to east-facing leaves increased more quickly 
290than did sap flow rates to west-facing leaves. East-facing leaves had sap flow rates of 
2911 .5 cm hr 1 on average 26 minutes before sap flow rates to west-facing leaves reached 
292the same threshold (t = 5.67, df = 23, P < 0.001 ; Figure 1). In addition, on 18 out of 25 
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293days, west-facing leaves reached their daily peak sap flow rate later in the day than 
294east-facing leaves. However, sap flow rates to east-facing leaves did not decline any 
295earlier in the evening than west-facing leaves, and east-facing leaves generally had 
296higher nighttime sap flow rates than west-facing leaves, perhaps indicative of greater 
297refilling. 

298 Patterns of sap flow to leaves and leaflets also differed. Leaflets generally had 
299lower sap flow rates than leaves, and the daily integrated sap flow through petioles and 
300petiolules differed significantly (t = 7.42, df = 24, P < 0.001 ; Figure 2). While water was 
301withheld for one week, daily maximum sap velocities for both leaves and leaflets 
302declined such that leaflet sap flow rates were about half of those to leaves (Figure 2a). 
303On the day immediately following re-watering, leaves and leaflets had almost equivalent 
304sap flow rates, which continued to increase on subsequent days despite declining daily 
305maximum VPD during these days. Daytime integrated sap flow to leaves and leaflets 
306was negatively correlated with mean VPD (Figure 3), both when including all days and 
307when the last five days of the drought treatment were excluded (Table 1). There was a 
308significant, negative relationship only between nighttime integrated sap flow of leaflets 
309and mean nighttime VPD, but only when all data, including the drought days, were 
310included. There was no relationship between nighttime VPD and integrated sap flow for 
31 Heaves. Maximum v h for leaves occurred at a slightly higher VPD than it did for leaflets 
312(2.21 kPa vs. 2.06 kPa; grey symbols in Figure 3). 

313 Patterns of sap flow hysteresis can be grouped into two classes, exemplified by 
314data from two days from the same leaf shown in Figure 4. Data in Figure 4 are 
315consistent with the patterns seen for other sensors on other days. The first type of 
316hysteresis pattern is denoted by clockwise hysteresis in the relationship between v h and 
317VPD (Figure 4a). On this day, hysteresis in the v h -PAR relationship was also clockwise 
318(Figure 4b). The second type of hysteresis is defined by an intersected loop, or figure- 
319eight pattern, in the relationship between v h and VPD (Figure 4d). On this day, the v h - 
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320PAR relationship had a counterclockwise pattern (Figure 4e). Nonetheless, the 
321relationship between PAR and VPD for these two days was similar, showing a 
322counterclockwise pattern for both days (Figure 4c,f). 

323 
324 

325DISCUSSION 

326 Leaf physiology responds rapidly to changes in the leaf microenvironment in 
327ways not previously appreciated (Zhang et al. 2013), and leaves may protect stems 
328from low water potentials that can lead to loss of xylem functioning (Sperry 1986, Hao et 
329al. 2008, Chen et al. 2009, 201 0, Johnson et al. 201 1 , Bucci et al. 201 2, Zhang et al. 
3302013). As a result, there has been burgeoning interest in the diurnal variability of leaf 
331hydraulic functioning (Brodribb and Holbrook 2004, 2007, Johnson et al. 2009, Johnson 
332et al. 2011 , Wheeler et al. 2013). The new techniques for measuring leaf-level sap flow 
333used here may become critical in quantifying diurnal variability in leaf functioning. Our 
334results highlight how water use can differ significantly even between adjacent leaves, 
335further justifying the need for fine-scale measurements like ours in quantifying leaf 
336responses to environmental drivers. Notably, our sap flow measurements on individual 
337leaves suggest that there may be important functional linkages between hydraulic 
338architecture and water use. While a number of studies have shown these linkages for 
339stems, substantially fewer attempts have been made to connect leaf hydraulic 
340architecture to water use under natural conditions. 

341 East- and west-facing leaves showed a number of differences in their patterns of 
342daily sap flow (Figure 1a). Midday depression of v h occurred in leaves on both sides of 
343the plant, although not always at the same time of day, probably due to a combination of 
344factors including differences in leaf energy balance and differences in the temporal and 
345spatial dynamics of leaf water potential changes. Furthermore, aspect significantly 
346affected the timing and rates of sap flow. Sap flow rates to east-facing leaves 
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347increased, on average, 26 minutes earlier in the morning than they did to west-facing 
348leaves, due to earlier increases in leaf-level PAR to east-facing leaves than to west- 
349facing leaves. Despite these differences in timing, leaf hydraulic conductance of east- 
350and west-facing leaves may be similar if higher transpiration rates in east-facing leaves 
35 tare accompanied by greater declines in leaf water potential. Leaf microclimatic 
352conditions can cause substantial differences between even adjacent leaves on the 
353same branch (Figure 1), which can influence the dynamics of branch sap flow (Burgess 
354and Dawson 2008). How much influence water use by one leaf may have on water use 
355by another, adjacent leaf is likely related to xylem hydraulic architecture. In addition to 
356flowing longitudinally from roots to leaves, water may also flow laterally within a stem 
357(MacKay and Weatherley 1973, James et al. 2003, Schulte and Costa 2010). The 
358degree of this lateral flow varies among species and results from lateral connections 
359between adjacent xylem vessels. Highly sectored xylem leads to close coupling of 
360water uptake by roots on one side of the plant and water use by leaves on the same 
361side of the plant. In this case, adjacent leaves on different sides of the plants would 
362draw upon largely different pools of water in the stem. High sectoriality in xylem 
363architecture allows plant parts to function independently, such that branches or leaves 
364may compete little for water (Brooks et al. 2003, Orians et al. 2005). In contrast, highly 
365integrated xylem (low sectoriality) leads to tighter hydraulic linkages between adjacent 
366leaves on different sides of the stem axis. While we do not know how well integrated 
367the xylem of adjacent leaves in T. rosea may be, orthostichous leaves (vertically aligned 
368along the shoot axis) generally have more interconnected vasculature than do non- 
369orthostichous leaves (those on different sides of a shoot; Watson and Casper 1984, 
370Orians et al. 2005). Thus, adjacent east- and west-facing T. rosea leaves probably 
371function more independently than would two east-facing, orthostichous leaves. Sap 
372flow to individual leaves varies among leaves on the same branch, and the magnitude 
373of this variation may itself vary among species depending on xylem architecture. 
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374 Patterns of sap flow through petioles were similar to patterns observed for 
375petioles of other tropical species (Roddy and Dawson 2012, 2013), but, in some cases, 
376different from patterns observed for main stems. On some days, patterns of hysteresis 
377in the i/„-VPD relationship were similar to those seen for main stems of canopy trees 
378(Meinzer et al. 1997, O'Grady et al. 1999, 2008, Zeppel et al. 2004; Figure 4a). In this 
379type of hysteresis, v h was higher in the morning than in the afternoon for a given VPD, 
380creating a clockwise loop in the relationship between v h and VPD. On days when the w 
381VPD relationship showed a clockwise loop, the v h -PAR relationship also had a clockwise 
382hysteresis loop (Figure 4b). In contrast, for main stems, a clockwise wVPD loop is 
383normally accompanied by a counterclockwise wPAR loop (Zeppel et al. 2004). On 
384days with this first type of hysteresis, v h was higher in the morning than in the afternoon, 
385with maximum daily v h occurring closer in time to peak PAR than to peak VPD. On 
386these days, PAR peaked early in the day, saturating stomatal conductance and leading 
387to high v h even when VPD was moderate. However, we observed a second type of 
388hysteresis, characterized by a counterclockwise loop in the v h -PAB relationship (Figure 
3894e) that, unlike for main stems, was accompanied by a markedly different relationship 
390between v h and VPD. When the v h -PAR relationship exhibited counterclockwise 
391 hysteresis, the wVPD relationship was characterized by an intersected loop, or figure- 
392eight (Figure 4d). Although this intersected loop has been reported previously for main 
393stems (Meinzer et al. 1999, O'Grady et al. 1999, 2008), its meaning has not been fully 
394discussed or understood. This pattern occurred when afternoon v h was higher than 
395morning v h , causing maximum daily v h to occur closer in time to peak daily VPD than to 
396peak daily PAR. If morning transpiration is low and does not result in substantial water 
397potential declines, then v h may peak in the afternoon when VPD peaks, as occurred on 
398the second day shown. Why morning sap flow on this day was so low remains unclear 
399but may be due to low water potentials, which we did not measure. Regardless, this 
400second type of hysteresis exhibiting an intersected loop requires (1) a bimodal peak in 
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401the daily v h pattern (i.e. midday depression of gas exchange, which commonly occurs in 
402tropical species) and (2) maximum daily v h to occur in the afternoon. Because slight 
403midday depression of v h occurred on both days shown in Figure 4, midday depression 
404alone may not lead to the intersected loop hysteresis. The second day did, however, 
405have a higher afternoon VPD both in absolute terms (maximum of 4.77 compared to 
4063.44 kPa) and relative to PAR (Figure 4c,f), which was probably partly responsible for 
407increased afternoon transpiration. 

408 The most probable cause underlying such different patterns of hysteresis for 
409individual leaves and for main stems is likely to be a matter of scale. Sap flow through 
410stems integrates the individual sap flow responses of many leaves in drastically different 
411microclimat.es. One of the most obvious sources of within-canopy variation is between 
412different parts of a plant canopy that undergo different diurnal patterns of incident PAR, 
413yet measurements on main stems ignore most of this within-canopy variation. In the 
414present study, leaf aspect influenced patterns of sap flow, and previous studies on 
415branches have shown that aspect influences both absolute rates of sap flow and the 
416timing of peak sap flow within the day (Steinberg et al. 1990, Akilan et al. 1994, Martin 
417et al. 2001, Alarcon et al. 2003, Burgess and Dawson 2008; Figure 1). Time lags 
418between daily peaks of sap flow for east- and west-facing branches of large trees would 
419result in different patterns of hysteresis depending on branch aspect (Burgess and 
420Dawson 2008), and these patterns for branches may be similar to the second type of 
421hysteresis (the figure-eight) we report for individual leaves. By measuring incident PAR 
422to each leaf, we attempted in the present study to account for some of the variation in 
423leaf microclimate that influences sap flow. However, we still ignored some important 
424factors, such as leaf temperature and its effects on leaf saturation vapor pressure and 
425the vapor pressure gradient (VPG) driving transpiration. This may be an acceptable 
426oversight because atmospheric humidity has a greater impact on stomatal conductance 
427than does leaf temperature (Fredeen and Sage 1999, Mott and Peak 2010). In addition 
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428to microclimatic variation, leaves and stems differ in their hydraulic architecture, which 
429could influence sap flow patterns and hysteresis. Leaf water balance changes rapidly 
430as efflux and influx of water vary asynchronously on the timescale of seconds (Sheriff 
43 land Sinclair 1973, Sheriff 1974). Water balance of stems may not change as rapidly, 
432however, because of the compensatory effects of having numerous parallel pathways 
433for water entry and loss. Thus, transpiration and sap flow may vary over much shorter 
434timescales for leaves than for stems. Examining and quantifying sap flow hysteresis 
435may provide new insights into hydraulic functioning in response to various abiotic 
436factors influencing transpiration (Zeppel et al. 2004, Pfautsch and Adams 2013, Roddy 
437and Dawson 2013). 

438 Integrated daily plant water use, as measured by sap flow, generally increases 
439with increasing mean and maximum daily VPD for plants from a wide variety of habitats, 
440including canopy trees and shrubs (e.g. Zeppel et al. 2004, Pfautsch and Adams 2013, 
441Skelton et al. 2013). However, in our experiment integrated daily leaf water use 
442decreased with increasing mean daily VPD (Figure 3), whether days of declining soil 
443water content were included in the analysis or not (Table 1). There was a significant 
444negative relationship between integrated nocturnal sap flow and VPD for leaflets, but 
445not for leaves, although this relationship was driven by very low nighttime sap flow 
446during the drought (Table 1). For both leaves and leaflets, the VPD at which maximum 
447daily v h occurred was remarkably well conserved across days (2.21 and 2.06 kPa, 
448respectively) and was, interestingly, the same whether maximum v h occurred in the 
449morning or in the afternoon (Figure 4a,d). These patterns opposite to those seen in 
450main stems may result from higher than normal VPDs during our experiment. Leaves of 
451 T. rosea saplings may rarely encounter such high daytime VPD under natural 
452conditions, and stomatal sensitivity to VPD may be responsible for the negative 
453relationship we observed (Oren et al. 1999b). At VPD above ~2 kPa, instantaneous sap 
454flow rates often declined, consistent with stomatal closure to regulate transpiration rate 
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455and leaf water potential. For T. rosea saplings, the VPD at maximum v h was higher than 
456the VPD at maximum g s of other species, perhaps because of the higher than normal 
457VPD during our experiment and the time lag between reaching maximum g s and 
458maximum v h due to hydraulic resistance. 

459 In response to declining soil water availability, daily maximum v h declined for both 
460leaves and leaflets despite increasing VPD during this time. Rewatering caused an 
461immediate increase in leaflet v h , such that it almost equaled leaf v h (Fig. 2). 
462Nonetheless, integrated leaflet sap flow was, on average, 30% less than leaf sap flow. 
463Assuming the ratio of leaflet area to conduit cross-sectional area (LA:SA ratio) is the 
464same for all leaflets, then instantaneous and integrated leaflet sap flow, as a fraction of 
465leaf sap flow, should be proportional to leaflet area. However, both instantaneous and 
466integrated leaflet sap flow were higher than this prediction, probably because cross- 
467sectional conduit area of petiolules is lower than that of petioles. This could result from 
468a combination of conduit taper and differences in the number of conduits between ranks 
469(McCulloh et al. 2009, 2010). Although we did not measure conduit dimensions, our 
470results highlight the potential linkages between leaf hydraulic architecture and diurnal 
471 patterns of water use at different scales. As of yet, there has been remarkably little 
472effort to connect xylem structure-function relationships to continuous, sap flow 
473measurements of plant water use. 
474 
475 

476CONCLUSIONS 

477 Understanding leaf-level responses to abiotic conditions is critical for modeling 
478plant responses to future climate change. In the present study, we found that leaves 
479often exhibit sap flow responses to abiotic drivers that are notably different from 
480responses of stems, for two main reasons: (1) stems integrate over many leaves, each 
481 with their own microclimate and (2) moving the sap flow sensor closer to the sites of 
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482transpiration removes the confounding influence of capacitance distal to most stem or 
483branch sap flow sensors. Thus, sap flow measurements on main stems may not 
484accurately describe leaf-level processes. Furthermore, we found significant variation in 
485sap flow patterns between adjacent leaves that are related to differences in the leaf 
486microenvironment. Differences in sap flow between leaves and leaflets are likely due to 
487differences in hydraulic architecture that influence patterns of water use. Future 
488measurements of sap flow through petioles could better elucidate the biotic and abiotic 
489drivers of transpiration dynamics under natural growth conditions. 
490 
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661Figure legends 

662 

663Figure 1 . (a) Three days of sap flow for east- (dashed) and west-facing (solid) leaves. 
664Tick marks on the horizontal axis indicate at midnight, (b) Boxplot of the time lag in sap 
665flow rates between east- and west-facing leaves. 

666 

667Figure 2. (a) Eight days of sap flow for leaves (solid) and leaflets (dashed). The first five 
668days corresponded to the end of a week without watering. Tick marks on the horizontal 
669axis indicate at midnight, (b) Boxplot of the daily time-integrated sap flow through 
670leaves and leaflets. 

671 

672Figure 3. The relationship between daily integrated sap flow and mean daily VPD for 
673leaves (triangles) and leaflets (circles). Black points represent the last five days during 
674the drought treatment. Grey points at the bottom mark the mean VPD (and standard 
675error) at which maximum daily sap flow occurred for leaves and leaflets. 

676 

677Figure 4. The pairwise relationships between v h , VPD, and PAR for two days (top and 
678bottom rows). (a,d) The v„-VPD relationship differed between the two days, as did the 
679(b,e) i/ h -PAR relationship. (c,f) However, the VPD-PAR relationship was approximately 
680the same for the two days. 
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681 

682Table 1 . Summary statistics for the linear regressions between integrated sap distance 
683and mean VPD for leaves and leaflets in the day and in the night, whether including 
684data from the week of drought or not. 

685 
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Figure legends 

Figure 1 . (a) Three days of sap flow for east- (dashed) and west-facing (solid) leaves. 
Tick marks on the horizontal axis indicate at midnight, (b) Boxplot of the time lag in sap 
flow rates between east- and west-facing leaves. 

Figure 2. (a) Eight days of sap flow for leaves (solid) and leaflets (dashed). The first five 
days corresponded to the end of a week without watering. Tick marks on the horizontal 
axis indicate at midnight, (b) Boxplot of the daily time-integrated sap flow through 
leaves and leaflets. 

Figure 3. The relationship between daily integrated sap flow and mean daily VPD for 
leaves (triangles) and leaflets (circles). Black points represent the last five days during 
the drought treatment. Grey points at the bottom mark the mean VPD (and standard er- 
ror) at which maximum daily sap flow occurred for leaves and leaflets. 

Figure 4. The pairwise relationships between v h , VPD, and PAR for two days (top and 
bottom rows). (a,d) The wVPD relationship differed between the two days, as did the 
(b,e) V/j-PAR relationship. (c,f) However, the VPD-PAR relationship was approximately 
the same for the two days. 

Table 1 . Summary statistics for the linear regressions between integrated sap distance 
and mean VPD for leaves and leaflets in the day and in the night, whether including 
data from the week of drought or not. 
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